Fits of the Thermodynamic Phase Plane

Abstract

Using the Weeks-Chandler-Andersen separation of the intermolecular potential we have fitted computer simulation data for Lennard-Jones system for the whole phase plane to the same form expression but in two different ways: locally and globally. We compare the efficacy and the exactness of both methods.     

Molecular dynamics simulation study of probe diffusion in liquid n-alkanes

We performed molecular dynamics simulations for the probe diffusion and friction dynamics of Lennard-Jones (LJ) particles modelled for methyl yellow (MY) in liquid n-alkanes of C₁₂–C₂₀₀ at temperatures of 318, 418, 518 and 618 K. Two LJ particles are chosen: MY1 with a mass of 114 g/mol, LJ parameters of σ = 4.0 Å and ? = 0.4 kJ/mol, and MY2 with a mass of 225 g/mol, σ = 6.0 Å and ? = 0.6 kJ/mol. We observed a clear transition in the power law dependence of MY2 diffusion on the molecular weight of n-alkanes at lower temperatures of 318 and 418 K. The sharp transitions occur near n-dotriacontane (C₃₂). However, no such transition is found for MY1 at all the temperatures and for MY2 at higher temperatures of 518 and 618 K. We also calculated the friction constants of both MY probe molecules in liquid n-alkanes. For the larger probe molecule (MY2), at lower temperatures, a large deviation of slope from the linear dependence of the friction of MY2 on the chain length of n-alkane is observed, which indicates a large reduction of friction in longer chains when compared with the shorter chains, enhancing the diffusion of the probe molecules (MY2).     

Study on the hydrolysis mechanism of phosphodiesterase 4 using molecular dynamics simulations

We carried out NPT molecular dynamics simulations in an explicit solvent to better understand the mechanism of cyclic adenosine monophosphates (cAMP) hydrolysis by phosphodiesterase 4 (PDE4) enzyme on atomic details and to obtain information on the dynamics characteristic of the catalytic domains of PDE4. In analyzing the water hydrogen-bond network around the active site, we also showed the importance of water in drug–protein interactions. In addition, we reported the characteristics of the hydration pattern and the dynamic distance distribution around the interesting residues. The results indicated that Asp318 plays the role of a general base that can activate water molecule for nucleophilic attack on cAMP. As expected, His160 plays the role of a proton donor for cAMP.     

How good are comparative models in the understanding of protein dynamics?

The 3D structure of a protein is essential to understand protein dynamics. If experimentally determined structure is unavailable, comparative models could be used to infer dynamics. However, the effectiveness of comparative models, compared to experimental structures, in inferring dynamics is not clear. To address this, we compared dynamics features of ~800 comparative models with their crystal structures using normal mode analysis. Average similarity in magnitude, direction, and correlation of residue motions is >0.8 (where value 1 is identical) indicating that the dynamics of models and crystal structures are highly similar. Accuracy of 3D structure and dynamics is significantly higher for models built on multiple and/or high sequence identity templates (>40%). Three-dimensional (3D) structure and residue fluctuations of models are closer to that of crystal structures than to templates (TM score 0.9 vs 0.7 and square inner product 0.92 vs 0.88). Furthermore, long-range molecular dynamics simulations on comparative models of RNase 1 and Angiogenin showed significant differences in the conformational sampling of conserved active-site residues that characterize differences in their activity levels. Similar analyses on two EGFR kinase variant models highlight the effect of mutations on the functional state-specific αC helix motions and these results corroborate with the previous experimental observations. Thus, our study adds confidence to the use of comparative models in understanding protein dynamics.     

Crystallographic structure versus homology model: a case study of molecular dynamics simulation of human and zebrafish histone deacetylase 10

Abstract

Histone deacetylase (HDAC) 10 has been implicated in the pathology of various cancers and neurodegenerative disorders, making the discovery of novel inhibitors of the isoform an important endeavor. However, the unavailability of crystallographic structure of human HDAC10 (hHDAC10) hinders structure-based drug design effort. Previously, we reported the homology modeled structure of human HDAC10 built using the crystallographic structure of Danio rerio (zebrafish) HDAC10 (zHDAC10) (Protein Data Bank (PDB) ID; 5TD7, released on 24 May 2017) as a template. Here, in continuation with our study, both hHDAC10 and zHDAC10, and their respective complexes with trichostatin A (TSA), quisinostat, and the native ligand (in 5TD7), 7-[(3-aminopropyl)amino]-1,1,1-trifluoroheptane-2,2-diol (PDB ID; FKS) were submitted to 100?ns-long unrestrained molecular dynamics (MD) simulations. Comparative analyses of the MD trajectories revealed that zHDAC10 and its complexes displayed higher stability than hHDAC10 and its corresponding complexes over time. Nonetheless, docking of active and inactive set molecules revealed that more reliable conformations of hHDAC10 could be obtained at an extended time period. This study may shed more light on the reliability of hHDAC10 modeled structure for use in selective inhibitor design.

Communicated by Ramaswamy H. Sarma     

Developing and testing alien species indicators for Europe

Alien species indicators provide vital information to the biodiversity policy sector on the status-quo and trends of biological invasions and on the efficacy of response measures. Applicable at different geographical scales and organizational levels, alien species indicators struggle with data availability and quality. Based on policy needs and previous work on the global scale, we here present a set of six alien species indicators for Europe, which capture complementary facets of biological invasions in Europe: (a) an combined index of invasion trends, (b) an indicator on pathways of invasions, (c) the Red List Index of Invasive Alien Species (IAS), (d) an indicator of IAS impacts on ecosystem services, (e) trends in incidence of livestock diseases and (f) an indicator on costs for alien species management and research. Each of these indicators has its particular strengths and shortcomings, but combined they allow for a nuanced understanding of the status and trends of biological invasions in Europe. We found that the scale and impact of biological invasions are steadily increasing across all impact indicators, although societal response in recent years has increased. The Red List Index is fit-for-purpose and demonstrates that overall extinction risks (here shown for amphibians in Europe) are increasing. Introduction pathway dynamics have changed, with some pathways decreasing in relevance (e.g., biological control agents) and others increasing (e.g., horticultural trade) providing a leverage for targeted policy and stakeholder response. The IAS indicators presented here for the first time on a continental basis serve as a starting point for future improvements, and as a basis for monitoring the efficacy of the recent EU legislation of IAS. This will need a better workflow for data collection and management. To achieve this, all main actors must work toward improving the interoperability among existing databases and between data holders.     

Identification of dengue viral RNA-dependent RNA polymerase inhibitor using computational fragment-based approaches and molecular dynamics study

Dengue is a major public health concern in tropical and subtropical countries of the world. There are no specific drugs available to treat dengue. Even though several candidates targeted both viral and host proteins to overcome dengue infection, they have not yet entered into the later stages of clinical trials. In order to design a drug for dengue fever, newly emerged fragment-based drug designing technique was applied. RNA-dependent RNA polymerase, which is essential for dengue viral replication is chosen as a drug target for dengue drug discovery. A cascade of methods, fragment screening, fragment growing, and fragment linking revealed the compound [2-(4-carbamoylpiperidin-1-yl)-2-oxoethyl]8-(1,3-benzothiazol-2-yl)naphthalene-1-carboxylate as a potent dengue viral polymerase inhibitor. Both strain energy and binding free energy calculations predicted that this could be a better inhibitor than the existing ones. Molecular dynamics simulation studies showed that the dengue polymerase–lead complex is stable and their interactions are consistent throughout the simulation. The hydrogen-bonded interactions formed by the residues Arg792, Thr794, Ser796, and Asn405 are the primary contributors for the stability and the rigidity of the polymerase–lead complex. This might keep the polymerase in closed conformation and thus inhibits viral replication. Hence, this might be a promising lead molecule for dengue drug designing. Further optimization of this lead molecule would result in a potent drug for dengue.     

Hard-sphere/pseudo-particle modelling (HS-PPM) for efficient and scalable molecular simulation of dilute gaseous flow and transport

Continuum methods are not accurate enough for flows at high Knudsen numbers, whereas rigorous molecular dynamics (MD) methods are too costly for simulations at practical dimensions. Hard-sphere (HS) model is a simplified MD method efficient for dilute gaseous flow but is of poor parallelism due to its event-driven nature, which sets a strong limitation to its large-scale applications. In this work, pseudo-particle modelling, a time-driven modelling approach is coupled with HS model to construct a scalable parallel method capable of simulating flows and transport processes at high Knudsen numbers without losing necessary molecular details in describing their macro-scale behaviours. The method is validated in several classical simulation cases and its performance is evaluated to be favourable. To demonstrate the potential applications of this method, we also simulate the diffusion of small molecules in multi-scale porous media which is related to catalysis, material preparation and micro chemical engineering in the long term.     

A model of litter harvesting by the Western Australian wheatbelt termite,Drepanotermes tamminensis (Hill), with particular reference to nutrient dynamics

A series of papers have been published which describe the influence of vegetation and soil type on the Western Australian wheatbelt termite,Drepanotermes tamminensis (Hill), and also on its litter harvesting levels and contribution to the soil nutrient budget. This paper integrates these findings by means of a computer simulation model. The model consists of three modules which respectively describe the dynamics of litter on the ground, the dynamics of litter within termite mounds and how these in turn influence nutrient loads within the habitat. The outputs of the model suggest that this litter harvesting termite plays an important role in the nutrient dynamics of the area and it provides an estimate of the unmeasured variable, litter consumed in mounds by termites, which is consistent with measurements for other termite species with similar feeding habits.